#include "slepdef.h"
#include "testDefs.h"
#include "slep.h"
#include <stdio.h>
#include <ctype.h>
#include "testblas1.h"

void exampleLogisticR();
void TestSlep_daxpby();
void exampleNNLogisticR();

int main1(){
	//TestAll();
	slep_matrix* mtrx;
	slep_vector* vctr;
	slep_matrix* Arnd;
	slep_vector* xrnd;
	slep_vector* yrnd;

	int i;
	//printf("EOF=%d, NULL=%d, \\n=%d, \\0=%d sizeofdouble=%d\n", EOF, NULL, '\n','\0',sizeof(double));

	//mtrx = slepReadDenseMatrix(SLEPWD "A.mtrx");
	//vctr = slepReadDenseVector(SLEPWD "x.vctr");
	//slepVerboseWriteDenseMatrix(SLEPWD "B.mtrx", mtrx);
	//slepVerboseWriteDenseVector(SLEPWD "y.vctr", vctr);

	//Arnd = slepAllocDenseMatrix(20,10);
	//xrnd = slepAllocDenseVector(20);
	//yrnd = slepAllocDenseVector(10);
	//slep_randn(Arnd->data.mtrx,10*20,0.,1.);
	//slep_randn(xrnd->data.vctr,20,0.,1.);
	//slep_daAxpby(SlepTranspose,1.,Arnd,xrnd,0.,yrnd);

	//slepVerboseWriteDenseMatrix(SLEPWD "Arnd.mtrx", Arnd);
	//slepVerboseWriteDenseVector(SLEPWD "xrnd.vctr", xrnd);
	//slepVerboseWriteDenseVector(SLEPWD "yrnd.vctr", yrnd);

	//TestMEMMacros();
	//testRunLR1();
	//testdaAxpy();
	exampleLogisticR();
	exampleNNLogisticR();

	//TestSlep_daxpby();
	return 0;
}
